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1. Introduction 

The most important discovery, that everybody expects at the LHC, is the discovery of the Higgs boson. 
The exclusive process which has the best experimental signature, as far as we know, is diffractive Higgs 
production. Having two large rapidity gaps between the Higgs and the recoiled protons, this process has a 
minimal background from QCD processes without the Higgs. However, the total cross section for diffractive 
Higgs production turns out to be very small, namely about 3 fb (see detailed estimates by the Durham 
group, in Ref. [jl|). The largest cross section out of all the different processes for Higgs production, is 
the cross section for inclusive Higgs production, which reaches 40 - 50 pb (assuming that the Higgs mass 
M H » 100GeV/c 2 ) . 
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Figure 1: The associate W boson and Higgs produc- 
tion. 



Unfortunately, the background for this process 
turns out to be so large, that any experimental mea- 
surement of the Higgs in an inclusive process, looks 
rather problematic. In his talk at GGI WS (Florence 
2007), Hannes Jung made the suggestion of measur- 
ing the associate Higgs production, which is the in- 
clusive production of the Higgs together with the W 
boson. This process has a cross section which is about 
1.5% of that of the cross section for inclusive Higgs 
production, but this smallness could be compensated 
by its improved experimental signature. This process 
yields the production of the W boson and the Higgs, 
at more or less the same rapidity values (see Fig. |l|). 
The value for the cross section for this process, can 
be found only to within an accuracy of 30%, due to 



large uncertainties in the values of the structure functions. However, the ratio 



a{pp -> WH + X) 
a(pp -»■ W + X) 



3.5 x 10 



,-8 



(1.1) 



is known to within 7% accuracy. In this paper, we estimate the cross section of the main background 
process, namely, the associate production of the W boson and the bb pair, which has a mass close to Mff, 
and is produced at more or less the same value of rapidity. It is easy to see, that if we assume that the W 
boson and the bb pair are produced independently of the different parton showers, the ratio of the signal 
to the background can be estimated using the following equation 



R 



a{pp -> WH + X) 
a(pp -»■ W + [bb] +X) 



Ri 



X) 



This simple estimate for the cross section for bb production, gives 



a(pp - [bb] + X) oc 5^ 



and, using a to t = 110 mb (see Ref. [§]), we obtain 



M 



(1.2) 



H 



R ~ 9 -=-10 



;i.3) 



However, we see two sources for an increase of the bb production cross section. These are the energy 
growth of the cross section due to the gluon structure function, and the positive correlation between the 
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Figure 2: The exchange of the single BFKL Pomeron. 

W boson and the bb production, which has been seen at the Tevatron Q. The estimates of both sources, 
are the subject of this paper. 



This paper is organized as follows. In the next section, we calculate the cross sections for inclusive W 
boson and bb production, with a particular focus on the energy dependence of the cross section. 

The main section of this paper is the third, where we consider the correlation function for inclusive W 
boson and bb production. We show that this correlation function can be reduced to the calculation of the 
specific enhanced diagrams, for the BFKL Pomeron interaction. The detailed calculations of the enhanced 
diagrams are presented in this section, which include the integration over the entire kinematic region. In 
particular, the overlapping singularities are taken into account. Calculating this diagram, and comparing 
this calculation with the experimental data obtained at the Tevatron, we will be able to provide a reliable 



estimate in section 4 of the ratio of Eq. (1.2), at the LHC range of energies. In the conclusion we summarise 
our results. 



2. One parton shower contribution: one BFKL Pomeron exchange and inclusive cross 
section 



In this section, we will discuss the calculation of the inclusive cross section. Due to the factorisation theorem 
{|, H , and the AGK cutting rules || , only the production from one parton shower contributes to the single 
inclusive cross section. We begin with a discussion of the one BFKL Pomeron exchange amplitude. 

In the calculations that follow, let ki denote the momenta conjugate to the size of the interacting 
dipoles. qi are the momenta transferred along the Pomerons. The diagram of Fig. ^ shows the scattering 
of the dipole with rapidity Y, and transverse momentum k, off the dipole at rapidity Y = 0, which carries 
transverse momentum ko, due to the exchange of one Pomeron. In this notation, the amplitude of Fig. || 
has the expression m 
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+oo . 

A (IP) = Yl f d 19 {k,q; ln ) e^^ Y - Y '^h{ ln ) g(k ,q^ 
where h(j n ) = ( j n - -) A(7„) and A(7 n ) = (7 n (l-7n))~ 



where the conformal variable j n , is defined in Eq. Q2,2|) as 



(2.1) 



1 + n 

In = — ^ h ^ 
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Figure 3: Contour enclosing singularities for integra- 
tion over the conformal variable 7. 



In 



1 



7n 



n 



+ iv 



(2.2) 



The contour C consists of the imaginary 7 axis 
from ±ioo, and the semi circle at 00, to the left of 
the imaginary 7 axis, which encloses all singularities 
in the integrand of Eq. (|2.1[) (see Fig. ||) . Since 

9 (X<?; 7) 9 (4,g;7) « {k 2 /kiy for k 2 > k 2 , the 
integrand vanishes on the semi circle at ± ioo, such 
that it is sufficient just to replace 



c 



t +100 



(2.3) 



In evaluating the integral of Eq. (2T), and in 
particular for the first enhanced diagram, the calcu- 
lations that follow will be more economical if the integrals are expressed in terms of the variable v, instead 
of the variable 7 n , (where the relationship between 7„ and v is given in Eq. fl2.2|) ). For the integration 
limits e + ioo < 7 < e + % 00 (as e — > 0), then Eq. ( |2.2| ) gives the corresponding limits of integration for the 
variable i^as— 00 < 1/ < 00, and one should sum over all real positive integers n in Eq. fl2.2j). In this way 
in Eq. (E]l]), one can replace the integration over 7„ with the integration over v using the notation 



E 



00 

E 

n=— 00 



dv 



(2.4) 



The notation i^LA v i corresponds to the integration over the quantum numbers associated 

with the continuous unitary variable irreducible representations of SL ( 2 , C) , defined in Eq. ( |2.2| ). The 
energy levels are the SL ( 2 , C) eigenvalues of the BFKL kernel, given by 



a 5 {V(l) -» (^(7n))} 



a 5 {2V(l) -V(7n) -^(l"7n)} (2.5) 
a s (^(l)-i2e^fl^L±i + ti/H (2.6) 
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where ip (/) = d / df (lnT (/) ) , and g(^,<f;7nj, are the eigenfunctions. Thus, the one BFKL 
pomeron exchange amplitude shown in Fig. 0, has this expression in v notation 



A (LP) = J dugp (k, q;n = 0, z/) e ^(^,u)(Y-Y') y 2 = ^ ^ gp ^ q ;n = 0, -v^j 

X 1 dvgp (k,q;v) e^ Y - Y 'K 2 X(0,v) g P (k ,q;-v} (2.7) 



where in terms of v 

\(n,v) = ——— — — — (2.8) 

In Eq. ( |2.7| ), it is assumed that n = 0, in order to find the largest contribution at high energy. From 
now onwards, the notation co(v) = w(n = 0, v) will be used (see Eq. ( |2.6D ), and A(z^) = X(n = 0, v) 
and g(k,q;i^J = g (k, q; n = 0, vj . In Eq. (|2~7l), g P ( q, k; n = 0, v) = T P (k,q) g (q,k;n = 0,v) and 
g ( q, k; n = 0, v) has been calculated in Ref.||], which takes the form the form 



2 \ iv ] 

ot) f< -" ) (2 ' 9) 



, . /l 1.3 1. „ . \ ^ / 1 1.3 1., . (g*) 2 \ 

where F (v) = ^ - -,„, i - - W) 1 - ^ j ^ ^ - -«,, - - - w , 1 - j 

2 -2iu 2 T 2 (\-iv)T(iv) „ , s 2- 2 ^7r 2 r 2 (1 — iv) T (iv) 

and where C x (v) = — ^ / / '- C 2 (v) = =rr \ 2.10 

V ' -iv Y 2 (\+iv)T(-iv) y J -iv r 2 (1 + iv) T (—iv) y ' 

/,From this definition, one notes the following interesting property of C\ (v), namely 

d (v) d (-v) = ^ (2.H) 
v z 



In the limit that q <C k, then the function F(v) defined in Eq. ( [2.10 ) tends to unity, and Eq. (|2.9| ) 
reduces to 

g(k,q,n = 0,v) = ^1 (lek 2 )™ (2.12) 
Tp(k,q) is the so called impact factor, which is equal to 

T P (k,q) = E/ X\d 2 x t d 2 y t \^({x uyi })\ 2 e^-^ {l - e^~^) (2.13) 
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Figure 4: The inclusive production in the BFKL Pomeron approach. 

where ^ ({xi,y«}) is the wave function of the proton written in terms of the colourless dipoles, where 
Xi and yi are the coordinates of the dipole % (here \x{ — yi\ is the dipole size). It is convenient to calculate 
the amplitude for single pomeron exchange of Eq. ([2,71), for the case of q = 0, since we are considering 
here the contribution of the BFKL Pomeron to the total cross section. Hence, assuming that k^ <C k 2 , 
and using the expression of Eq. ( |2.9| ) for g ( q, k; n = 0, v) 



7T 4 



A (UP) = — T P (k, q = 0)T P (k , q = 0) P (k; k \Y - Y r ) 

(k 2 kl) 

w ^M(X-y) x(0,u) (2.14) 

It should be noted that P (k; ko\Y — Y') is the dimensionless amplitude for the single BFKL pomeron 
exchange, where as A ( IIP) is related to the cross section, and it has the dimensions of GeV -2 . The 
integration over v can be evaluated by expanding the BFKL function uj ( u) , around the saddle point v = 
(see Eq. (£§)) 

u{v) =u(0) - -i/V' (0) (2.15) 



Then the integration over v in Eq. (2.14), reduces to a Gaussian integral, and evaluating the v inte- 



gration, gives the expression for single Pomeron exchange as the expression 



[A ( If) Uo = ^ FP (fc ' 9 = Q) ^ q = 0) P(k;k \Y- Y') (2.16) 

(k 2 k 2 Y 

i 

where P (k; k \ Y - Y>) = ^ ( 0) y ( ^^r^yy ) " (2-17) 
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It should be stressed, that the value of v in this integral is small, namely (y = ln(£; 2 /&;Q)/(u;"(0)y) <C l), 
justifying our expansion in Eq. Q2.15 ). It should be noted, that in Eq. (2.14), contains as the initial condi- 



tion at Y — Y' = 0, the exchange of two gluons in the Born Approximation of perturbative QCD. Indeed, 
in this case, the contour in v can be closed on the singularities of A(n = 0, v), which are the poles of the 
second order at v = ii^. Considering the case when k 2 > /cq, we close the contour on the point v = +i i, 
and obtain the following expression 

P(k;k \Y-Y' -0) — )Mhi(k 2 /k 2 ) (2.18) 



which leads to 



vr 4 r P (k,q = 0) Tp (k Q ,q = 0) 



A (IIP) = ^^,*-^*>w,*-»j l n ( k 2 /k 2 ) (2.19) 



Eq. (pTlSj) is the answer for the Born approximation of perturbative QCD, (namely, for the case of 
two gluon exchange). For the case of calculating the inclusive cross section, we can use the kt factorisation 
theorem || , using the Mueller technique || , which allows one to reduce the calculation of the inclusive 
cross section to the case of single Pomeron exchange, using the optical theorem (see Fig. |], where the 
inclusive production of the W boson is shown in the Mueller technique ||). Indeed, according to the 
factorisation theorem, the inclusive production can be written as |T( 

da 



dY w 



J dxidx 2 J d 2 q 1 d 2 q 2 (j)(x 1 ,q 1 )(j)(x2,q2)cr(M w ,q 1 ,q2)5(x 1 X2S - Mw) (2.20) 



where s is the energy of the collision, and (j)(xi,qi) is the un-integrated structure function, which is 
related to the gluon structure function in the following way 



XiG( Xi ,q') = I dk 2 (P(xi,k 2 ) (2.21) 



where the Xi are equal to 



Xl =^e^; X2 = ^e"^; (2.22) 

where Yyy is the rapidity of the W boson in the center of mass frame (c.m.f.). The un-integrated 
structure function at high energy can be rewritten through the Pomeron exchange, while the hard cross 
section in Eq. ( p. 20 ) is calculated in the appendix. The integrals over q± and c/2 in Eq. ( 2.20 ), are convergent 



at q\ w My/ and q2 ~ M\y, and because of this, one can rewrite Eq. ( |2.20| ) in a simpler form, namely 
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da ind (ln(s/M fe y ;Y bb ,0 



dy, 



bb 



(2.23) 



a(M bb ) T P (k,q = 0)T P (ko,q = 0) J ^ P ( k,M bfi \Y -Y bb - ±6Y bb ) P ( ' k , M brb \Y bb - ±5Y, 



U2 h.2 
rv r\iQ 



6,6 1 bb 2 bb 



daind (In (f/Mgj ;Y W ;0) 
dY w 



(2.24) 



a(M w ) Tp(k,q = 0)T P (k ,q = 0) J P (k, M W \Y - Y w - -61%) P (ko, M W \Y W - \sY w 



In Eq. ( |2.24| ), a(M\y) is the squared amplitude of the quark hexagon contribution shown in Fig. ||, 
and a ( M bb ) is the squared amplitude of the quark subprocess contribution shown in Fig. [| Both have 
the dimensions of cross section, namely a{M\y) and a(M bb ) have dimensions of GeV -2 (see Eq. ( p. 27 ) 
and Eq. ( p. 29 ) respectively). Y = ln(s/m 2 ), and Y\y and Y bb are respectively the rapidity values in the 
laboratory frame of the W boson, and the center of mass of the quark anti-quark pair. SYyy = h^Mfy/m 2 ) 
and 5Y bb = ln(M 2 g/m 2 ) are the rapidity windows occupied by the W boson and the [bb] pair, respectively. 
The arguments in the single Pomeron amplitudes in Eq. ( 2.24 ), follow directly from Eq. ( |2.22| ). For P, we 
can use Eq. ( [2.14D , for example 



P[k,M w \Y-Y w --6Y w 



°° ±_ Mu) (Y-Y W ~ISY W ) ( M?> ^ "' 



,2iri 



1 w 



A: 2 



(2.25) 



Using Eq. (|2.18), one can see that at low energies, Y — > 5Yw, and Eq. (2.24) reduces to the expression 



~ (ln( ^ );1V;0) — T P (Jfe, q = 0) T P fa, q = 0) In (M 2 , A 2 ) In (M 2 ,/^ 2 ) o (M w ) (2.26) 

which is exactly the same as the expression for the cross section in the Born approximation of pertur- 
bative QCD, ( see Eq. fl2.2gj )), with k 2 = k$ = q 2 ± • , and Tp (k, q = 0) = as(k 2 ). It is easy to see, that 
Eq. ( 2.23| ) has the form of Eq. ( p. 30 ), in the low energy limit. 



The Mueller diagram in the Born approximation of pQCD for the inclusive W boson production, is 
shown in Fig. ^, and the formula for this cross section is given by Eq. ( [2,24 ). We need only to calculate 



a (Mw) and a(M bb ). The easiest way to derive these cross sections, is to calculate them in the Born 
approximation of perturbative QCD ( see Fig. |5| and Fig. ||). The calculations of the expressions for the 
cross sections for the processes of Fig. ||] and Fig. ||, are presented in appendices A-3 and A-4, and the final 
expressions for the amplitudes of these diagrams take the form 
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• for a (M w ,q 1 ,q 2 ) (see Fig. § and Eq. QA.l.lSQ ) 

a 2 (a 2 ) 

a(M w , qi ,q 2 ) = 2 SV y a{M w ) 



Recall that the Fermi coupling Gf has the dimensions of GeV -2 , so that inspection of Eq. Q2.27D , 
shows that a (M\y) has the dimensions of the cross section. The factor in front of a(Mw), (in 
the first line of Eq. ( |2.27| )), is taken into account in the BFKL Pomeron exchange in Eq. fl2,24| ). 
Integrating over the transverse momenta q\ t ± and q 2j ± on the RHS of Eq. ( |2.27| ), one obtains the 
following expression for the Born amplitude of Fig. [|. 



cr(M w ,qi,q 2 ) = a|(g 2 )7r 2 ln 2 



<r(M w ) 



(2.28) 



for a gi, ga) (see Fig. | and Eq. flOp ) 
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a 2 (a 2 ) 

a(M bh q u q 2 ) = 2 5W / a (M bb ) 

*1,± 12,1. 

(N 2 -l\ 2 la 2 1 (M 2 \ , . 

where a (M hb ) = — c - £ In — ^ 2.29 

Integrating over the transverse momenta q\± and (7*2, l on the RHS of Eq. (|2.29|) , one obtains the 
following expression for the Born amplitude of Fig. @. 



a(M bb ,q 1: q 2 ) = a 2 (q V In 2 U (M bb ) (2.30) 

\?l,min / 

3. Two parton shower contribution to the background for W boson and Higgs produc- 
tion 

3.1 The simplest diagram 

For the associate W boson and Higgs production, we expect a small background, since the main process 
which creates such a background, namely, 

p + p^bb + W + X (3.1.1) 

leads to a negligible contribution if both the W and the bb pair are produced from one parton shower, at 
Y bb « Yw- Here, Y bb is the value for the rapidity of the quark anti-quark pair, namely, Y bb = ^(Y b + Y b ) and 
Yyy is the Higgs boson rapidity. Indeed, in one parton shower, the typical difference (Ay) in the rapidity 
between the two emitted partons, are larger (or equal to) l/as 1. Therefore, keeping \Y bb — Yw\ < l/ass 
we have a suppression for the production of the W boson and the quark anti-quark pair. 

However, there exists a process for the production of the quark anti-quark pair and W boson production, 
from two different parton showers (see Fig. [?]). In this process, we do not expect any suppression of the 
quark anti-quark pair, and W boson production, at the same values of rapidity. However, this has its own 
suppression, which is related to the small probability associated with having two parton shower processes. 
However, at high energies, this suppression is not actually very strong, as we will demonstrate below. 

The general expression for the diagram of Fig. [7| has the following form 



A(Fig.§ 



d 2 q 
(2n) 2 



V 2 (pp -» JPJP; q) 



da incl (ln(s/M 2 -);Y bb ; -q 



dY 



bb 



doing (ln(g/M^ 
dY w 



(3.1.2) 
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V (pp — ► PIP) in Eq. (|3.1.2| ) and Eq. (|3.1.3|) below is the vertex of the interaction of the two Pomerons 
with the proton (see Fig. 0). In Eq. ( 3.1.2j ), doi nc i/dy is the cross section for the inclusive production 
of the W boson, and/or the quark anti-quark pair, but not at q = 0. However, at high energy and for 
hard processes, we can assume that the value of the typical q is small, and it is of the order of the typical 
momenta inside of the proton. At high energy, the mechanism for inclusive production is related to the 
emission of one parton shower, which can be described by the DGLAP evolution [11 1, or by the BFKL 
Pomeron exchange 0. For both cases, it has been proved [12, [D], || that we can safely put q = in our 
evaluation of the integral of Eq. (3.1.2). Therefore, 



dY b i Y w J (2-kY 



d(T r 



ind. (Hs/Ml- b );Y b - b ;0) dcT mcl (Hs/M^);Yw;0) 



2a, 



eff 



(3.1.3) 



Fortunately for us, the value of a e ff has been measured by the CDF collaboration at the Tevatron ||, 
with the result <J e ff = 14.5± 1.7±2.3 mb. Using the value of the predicted cross section for bb production, 



taking M bi = M H = 100 GeV/c 2 @ for Y bi = 0, and at the LHC energy, namely, 
2 x 10 _6 mb*, one can obtain from Eq. ( 3.1.3[ ) (comparing it with Eq. (1.1)), that 

d<Jbackground(Eq. (|3.1.3|) ) / dy h l\Y ll= 



da mcl (ln(s/Ml- b );Y b - b ;0) 



da H production/ dY}{ \ Y H =0 



2 in the LHC energy range (3.1.4) 



*This value, together with a more detailed consideration of the value of the cross section for bb inclusive production, will 
be discussed in section 4 
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This ratio can be easily reduced using the characteristic property of the process of Fig. [7], namely that 
there are no correlations in the rapidity values between the W boson, and the quark anti-quark pair. In 
other words, the cross section for this process does not depend on the difference Yw — Y bb , while the cross 
section for the associate W Higgs production of Fig. |], has a maximum at Yw — Y bb = 0. 

This result is encouraging, but we need to estimate the more complicated diagrams to derive the final 
conclusions (see Fig. |8|). 

3.2 The two parton shower production: enhanced diagram 

In the diagram of Fig. ||, the two parton showers are produced at the rapidity yi, and they merge back 
to one parton shower at the rapidity y 2 . Such a process can be reduced to the enhanced BFKL Pomeron 
diagram at high energy, as it is shown in Fig. ||. We would like to recall, that in this simplification, we use 
the unitarity constraint 



21m A 



el 



\A 



+ Gi r , 



(3.2.1) 



where A e i is the elastic scattering amplitude at energy s and impact parameter 6, while Gi n is the 
contribution of all inelastic processes. The scattering amplitude is purely imaginary at high energy, and it 
can be described by the exchange of one BFKL Pomeron. Neglecting \A e {\ in Eq. (3.2.1), which is possible 
for the range of energies asln(s / sq) <C ln(l/a|) , we obtain from Eq. ( 3.2.l| ) the relation 



2 A (IIP) = G m (s,b) (3.2.2) 

where A (IP) is the amplitude for one BFKL Pomeron exchange. Using this relation, we are able 
to replace the production of the parton showers with the exchange of the BFKL Pomerons. In this 
subsection, we will calculate the simplest enhanced diagram shown in Fig. ||, whose contribution is given 
by the expression 



Y 

I dYi 

'Yw+hSY w Jo 



2 

A(P^2P^P;Fig.§ = I h(i) d 7 f[<[ h( 7l ) drft I h(^) Hi <f ^(V) H 

dY 2 g P (k,q = Q;- 1 ) A -1 (7) e^ 7 ) ( y_Yl ) e^) (*-*W-§«V) 

x e ^(72) (Yi-Y bb -±5Y bh ) e w( 7 () (Y W -\8Y W -Y 2 ) ^(7^) (Y bb - \5Y bb -Y 2 ) ^(7') (Y 2 -0) 

xg P (ka,q = 0;^ A -1 (7') a(M w , -71, 7i) <r(M bl , -72,72) 

x Jd 2 q'T 31P (q = 0, -g"|7, 71, 72) T 3]P (q = 0, -q'\ - 7', -7^ -7Q 



(3.2.3) 



where h(^) and A(7) = (7(7 — 1)) have been defined in Eq. (2A). The contours C, C , C 



C2 , C[ and C 2 are similar to the contour shown in Fig. ^, which consists of the imaginary 7 axis from 
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Figure 8: The first enhanced Mueller diagram for two parton shower contribution to the background for the W 
boson and Higgs production 



±ioo, and the semi circle at oo to the left of the imaginary 7 axis, which encloses all singularities in the 



integrand of Eq. ( 3.2.3 ). In this case also, it is assumed that the integrand as a function of 7, 7', 71, 72, j[ 
and 7 2 , vanishes on the semi circle at 00, such that it is sufficient just to replace for each integration over 
the conformal variables 7/, in Eq. ( |3.2.3 ), 



C 



£ +200 



ill} = 71, 72,71) 72 



(3.2.4) 



Following the definition of Eq. ( p.2[ ), one can replace the integration limits for the 7; variables e — icc < 
7/ < e + ioo, with the corresponding limits of integration for the variables v\ {{v{\ = v, v' , u\, 1/2, v \i ^2), 
namely — cxd < v\ < 00. So in Eq. ( |3.2.3| ), one can change from the variable 7^ to the variable v\, and 
re-express the integration over 7^ as 



r p€ +200 p 

f d<yi -» / d7; = / 

J C J e— i 00 •/ it 



fia+oo 



(3.2.5) 



where a 



1/2 



e is such that the contour of integration is below all singularities in v\. Hence, 

du is understood to be a 



in v notation, Eq. ( 3.2.3 ) will be rewritten as follows, where the notation J 
shorthand for the contour integrals written in Eq. ( |3.2.5| ). 
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/OO _ f'OO f'OO f'OO p 

dvv 2 ]"[ / dvivfX^i) I d^vf \{y[) I dv'v' 2 J d 2 q 
-co i—\ 00 J —00 J — OO J 

Y - 

dY 2 g P (k,q = 0;-v) e <"(") O^i) e ^i){Y x -Y w -\SY w ) e ^u 2 )(Y 1 -Y b - b -^8Y b - b ) 



Y 

I dY 1 

'Y W + \SY W Jo 



Y bb ^ Y bb 



x ^K) {Y W -\&Y W -Y 2 ) e u>(v> 2 ) {Ybb -\5Y bb -Y 2 ) ^{v>) (Y 2 -0) gp ^ - = q. „A _^ ^ ^ 

x r 3]P {q = 0,-q'\u : u 1 ,u 2 ) r 3JP (q = 0, -q' \ - v', -v' x , -v' 2 ) (3.2.6) 



where a(M]y, — v\, v[) 

and a(M bb ,-v 2 ,v' 2 ) 



g(M w ,q = 0;-v 1 ) g(M w , q = 0; v[) a(M w ) 
g(M bl , q = 0; -v 2 ) g(M bb , q = 0; v' 2 ) a{M bb ) 



(3.2.7) 



where the expressions for the a {My/) and a (M bb ) are given in Eq. fE27p and Eq. ( ggg ), respectively. 
In Eq. ( p. 2. 6 ), we assume that IV > Y" fe g, and that q = 0, since we are considering the inclusive process 
given by Fig. ||, (see the left part of this figure). From Eq. ( 3.2.6j ), we see that it is possible to integrate 
over q 1 . This integration has been performed in Ref. ||, with the result 



, , 2-7TQ?. 

r 3 jp [q = 0,q \u,ux,u 2 ) = [ -j^- 



d k g(k,q = 0,v) g [k + -<f , q ', vi) g\k + -q", -q" , v 2 



1-/ 



/ 10 \ — 1/2+iv— iv\— IV2 



(1/2 + iv — iv\ — iv 2 ) 



+ function without singularities in v\ and v 2 

(3.2.8) 



In Eq. ( 3.2. 8; ), we restrict ourselves to the case of small v\ and v 2 ( v[ and v' 2 ), since as we will show 
below, the main contribution to the integral over rapidity values, stems from Y\ — Y 2 — » Y, and small values 



of vi and v 2 ( v[ and v' 2 ), contribute to the BFKL Pomeron exchange at high energy (see Eq. ( 2.15 ), and 



the discussion around this equation). From Eq. (|3.2.8|) , the result of integrating the product of the two 
triple pomeron vertices, which appear in the integrand on the RHS of Eq. ( |3.2.6 ), is 



J d 2 q'T 3IP (q = 0,q '\v ,v x ,v 2 ) T 3IP (q = 0,q'\ - v' , -v[ , -v' 2 ) = 



2ira 2 s 



( 1/2 + iv - ivi - iv 2 y ( 1/2 - iv' + iv[ + iv' 2 ) 



/\3 



x 5 ( v 



(3.2.9) 

v\ - v 2 - v + v[ + v 2 ) 



The Dirac delta function appearing on the RHS of Eq. ( 3.2.9 ), is absorbed in Eq. ( 3. 2. 6] ) by integrating 
over v' . In the course of this, the expression on the RHS of Eq. ( 3.2.9| ) reduces to 
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/ -2 \ 2 

( ^Wr) Ci(v) C l (v-v l -v 2 + v[ + v' 2 ) 2vr 9 ( 1/2 + iv - iv x - iv 2 y 3 ( 1/2 + iv - iv x - iv 2 )~ 3 . 
Hence, after evaluating the integration over v', and making use of the expression in Eq. ( fmi) , in the case 
where q = for the pomeron couplings, gives for the RHS of Eq. ( p. 2. 6 ) 



A(P -> 2JP -> iP;Fu/. 



27T a 



27ra 



2 \ 2 



T P (fc, g - 0) T P (fco, g = 0) a (M w ) a (M, 



WW 



/OO 2 /*00 /'OO 

d " 7/2 n / ^ 2 a ( ^) / d ^^ 2 a (^') K+^-^-^+f) 1 
-oo i = -^J — OO J —OO 

Y 



dY 2 g[k,q = 0;-v) e' 



wM (y-vi) p (u(^)(yi-ywr-i*y¥ir) ^(^(y-y^-iay^) 



x / cZYi 
'yvK+^yw 

x 5 (M W , -i*) <?(A%, i^) 3 (M 65 , -i/ 2 ) 5 (M b5 , i^) 

x Ci ( v) C\ ( i^i + v 2 — v[ — v 2 — v) (1/2 — iv + iz^i + iv 2 ) 3 (1/2 + iv — iv\ — iv 2 )~ 3 



x 9 (ko, q = Q;Vi + v' 2 — vi — v 2 + v 



(3.2.10) 



where from the expression of Eq. ( |212D , the product Ci C13 g can be recast in the form 



Ci (v) C\ ( -i/J - 1/2 + ^1 + f 2 - ^) g (k,q = 0; -z/J 5 [k ,q = 0; z/( + z/ 2 - v\ - v 2 + vj (3.2.11) 
— ► Ci ( v) C\ ( -i/) Ci ( i/( + z/ 2 - z/ a - z/ 2 + i>) Ci ( -i/J - i/ 2 + z> a + v 2 - v) 
{k 2 kl)- l,2 (lQk 2 y w (16 ^y^-i^-i^+iu 



Hence, it is obvious that using the relation of Eq. ( 2.11 ), the expression in Eq. ( |3.2.11 ) cancels the 
term v 2 ( v[ + v' 2 — v\ — v 2 + v) 2 in the numerator of the integrand of Eq. ( 3.2.10 ), to give 



A(P^2P^P; Fig. |) = ^ 8 S (k, M w , M bi , k ) 



x 



-OO 

Y 



i=l 

dY 1 



du W dviv 2 \(vi) / dv'^f \(v[) (16 k 2 ) lv 



2 \vu\+iv^~%u\—iv^ J rW 



Y - — - SY - 

" b 2 dYo e ^H( y -^i) e ^(u 1 )(Y 1 -Y w -l8Y w ) e c(u 2 )(Y 1 -Y b - b -^5Y b - b ) 



Iy w +±8y w Jo 

x e uj{u[) (Y W -\5Y W -Y 2 ) e u,K) (Y bb - \5Y bb -Y 2 ) e ^(u 1 +u 2 -u' 1 -u' 2 +u) (Y 2 -0) 

x g(M w , -v x ) g(M w , v[) g(M bh , -v 2 ) g(M bh , v' 2 ) ( 1/2 - iv + iv x + iv 2 y 3 ( 1/2 + iv - iv x - iv 2 )~ 



(3.2.12) 



where the factor 
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£ (*, M W , M 6 - 6! feo) = 2 vr 9 (^-^ J T P (k, q = 0)T P (k , q = 0) -j=^ a (M bb ) Mfa (M w ) 

(3.2.13) 

Substituting for g appearing in Eq. ( |3.2.12j ) for the form given in Eq. ( |2.9D , the expression on the RHS 
of Eq. ( |3.2.120 becomes 



A{P^2P^IP ] Fig. D = 7T 16 £ (fc, Af»r, M bb , k ) 

2 



X 



/oo /*oo z 1 oo 

d^n/ d^A(^) / d^A(^) (lGfc 2 )^ ( 16A .2)-i+- 2 -M-^+- 
-OO .• i J —oo J —oo 



1=1 

,-Y 

X 



x e c,(^)(y w -i5y w -y 2 ) e ^K)(y bE -i5y i)6 -y 2 ) e a;K+^- l , 1 -^+ ! ,)(y 2 -o) (i6M^)"^~ W1 (l6M 2 -) iiV ~ ii/2 

x (l/2-ii/ + ^i + ii/ 2 )~ 3 (1/2 + izz-M/! -^ 2 ) -3 (3.2.14) 



Note that due to the Vi and u- (i = 1 , 2) in the denominator of Eq. ( |2.g| ) , inserting this expression for 
the g's into Eq. ( |3.2.12| ), the vf and v'^ in the numerator of Eq. ( |3.2.12[ ) reduce to just factors of Vi and v[ 
in the numerator of Eq. ( |3.2.12 ). Since u(ui) > contains a minimum at V[ = 0, the integration over Y\ 



and Y2 is actually convergent, either at Y — Y\ < l/u;(0) ( Y% < /oj(0)) or at Y\ — Yw — \5Yw < 1/cj(0) ( 
Y bb - \5Y bb -Y 2 < 1/^(0)). The maximal increase with energy stems from the region where Y\ — > Y and 
Y% — ► 0. Indeed, in this region the integrand in Eq. ( |3.2.12j ) is proportional to 

e u{ux){Y-Y w -\5Y w ) e u{u. 2 ){Y-Y b - b -\5Y b - b ) x e ^(u[) (Y W -\5Y W ) ^(u' 2 ) (Y bb -\5Y bb ) (3.2.15) 

Using Eq. ( ^3.2. 15| ) , and evaluating the integrations over all the z/'s in the saddle point approximation, 
we see that 

A(P^21P^1P; Fig. §) oc exp L(0) Uy - )ff w - ^5Y bb ^j } (3.2.16) 
where the typical values of v \ are 



^ HM^/kj) ^ HMl/kp 

1 11 ~ 2(y - Y W - \5Y W ) ' 1 21 ~ 2(Y - Y bb - \5Y bb ) 

1 11 20V - * 1j N 2(y 65 - ^) « h (3 ' 2 - 17) 
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This can be seen more clearly in the argument that follows, where considering the region of integration 
close to the point v — ► 0, one obtains the delta function 8 [Y — Y\ + Y2) in the rapidity variables. 

For small v\ and z/ 2 ( z/[ and u 2 ), we have three different kinematic regions of integration: (i) |iz/+l/2| <C 
v\ + ^2 + 1/2| < v[ + v' 2 ); (ii) \%v + 1/2| » z/i + z/ 2 flu/ + 1/2| » i/( + z/,); and (iii) |tV + 1/2| « z/i + f 2 
(|»V + 1/2| wi/ji j/ 2 ). 

3.2.1 Region (I): |u/ + 1/2| < i/i + i/ 2 (|u/' + 1/2| < i/£ + 

In the first kinematic region, we have I^jp — ► \/{iv\ + <if 2 ) 3 . Therefore, the integral is large for values of 
\v\ — > 1/2. Using Eq. ( |3.2.9D , the integrand in Eq. ( |3.2.6D could be simplified in this kinematic region in 
the following way. 



2na 
~N, 



2 n 2 



^reg.(i) (P^21P^ IP; Fig. § = 2vr 9 ( ) T P (k, q-0)T P (k ,q = 0) a (M w ) a (M a ) (3.2.1.1) 



x 



/oo 2 />oo /-oo />oo r/- _l 7, t/' \n2ii/\+2ii/2—2iu', —2iv' 



xg(k,q = 0;-v) / dY x dY 2 C^-n) e C"W (n-*V-£«V) Cn-V*-^) 

x g(M w , -m) g(Mw,v[) g(M a , -v 2 ) g(M bl , u' 2 ) g(k,q = 0; / 



The BFKL eigenvalues w(z/) and cj(z/), (see Eq. (El])) increase as — > —05/(1/2 ± viz/), and 
therefore, the integral over v takes the following form 



/ o- ex P -TJ^T- {Y-YM^-as — —2— = a s (Y - Y ± ) 3.2.1.2 

where t = 0:5/(1/2 ± iv) . The typical v in this integral, is about a$(Y — Yi), while we expect 
v\ ~ z/ 2 ~ 1/asY. Therefore, for region (i) |iz/ + 1/2 1 cannot be less than v\ + 1/2, which is in clear 
contradiction with our initial assumption. It means, that the first kinematic region does not contribute to 
the integral. 

3.2.2 Region (II): \iv + 1/2| > v x + v 2 {\iu' + 1/2| > u[ + v' 2 ) 



In the second kinematic region, \iu + 1/2 1 3> v\ + z/ 2 (|*z/ —I— 1/2 1 3> ^ + z/ 2 ) while z^i 
1/asY. In this region, Eq. ( ET2~3 ) reads 



^2 



J d 2 q'T 3] p(q = 0,q\u,ui,U2) T 3]P (q = 0,q'\ 



4) 



(3.2.2.1) 



N r . 



(l/2 + ii/) 3 (l/2-iz/) 3 



8{u - v') 



Taking into account Eq. ( |3.2.2.1[ ), the integration over v' is eliminated by the Dirac delta function 
5 ( v — v 1 ) , which appears in Eq. ( |3.2.2.1| ), such that inserting this into Eq. ( 3.2.6| ) gives the following 
expression 



A ieg . {ii) (lP^21P^ P;Fig.$) = (3.2.2.2) 

dv ( kV 



/oo j v /k \ * roo f-oo 

/ dY 1 \ 

Jy w +±sy w Jo 



i=l ' 

dYo e ^)i Y -Y 1+ Y 2 ) e w(u 1 )(Y 1 -Y w -^SY w ) ^(ya)(Xi-Y a -^6Y a ) 



x ay 1 ar 2 

Jy w +±sy w Jo 

x ^K) (Y W -\5Y W -Y 2 ) ^K) (Y bl -\5Y b - b -Y 2 ) ff(M ^ _^ ff(M ^ ^ g ^ _^ g{M ^ ^ 

If the following definitions are made, namely 

1 Myy 1 Mjt 

Yin = Y\ — Yy/ — -SYyy ; L V1 = „ ; Y U2 = Y\ — Y bb — -5Y bb ; L v% = —y ; 

1 Mr% 1 Af^ 

^ = Y W --SY W -Y 2] L v[ = -f; Y„, = Y a - -SY bb - Y 2 ; L v 



Then Eq. ( |3.2,2.2j ) can be written in the following, more compact form, 



(3.2.2.3) 



^reg.(u) (P^21P^ P; Fig. § = 

(I6vr 4 ) 4 £ (k, M w , M bi , k ) / dY x / dY 2 I v I Ul I U2 l' Ui l' V2 (3.2.2.4) 

Jy w +\sy w Jo 



,oo-i(e-l/2) e l^;0'-n+y 2 ) /fc2\^ 

where I U (Y-Y 1 -Y 2 ) = / <fo , 1/9W1/ . . , 3 r| (3-2.2.5) 

-oo-*(e-l/2) + V 2 ) (V 2 -^) / 



/oo 
du lVi e^ Y ^ (L u .) iUi (3.2.2.6) 
-oo 



i/i-fO d 



comparison with Eg. (|2.14|) =► J ((y L ) ^ — - — - — P (L u ; Y v . ) (3.2.2.7) 

am 

In Eq. ( |3.2.2.4| ), we evaluate first the integrals over and Yjj, which take the following form, where 
in region (ii), v\, u 2 , v[ and u' 2 are considered small. 



dY 1 / dY 2 I v I u I U2 1 I = / - - — — — -f 

W+f«V -/o y-00 («/ + i/2) 3 («/- 1/2) 3 V*; 2 / 

/°° r 00 7^0 L~ iv2 di/n r°° ■ 1 r°° vUL ? 2 dvL 

x ^(^)(K-y W '-i»Hr) e w(i<a)(y-n 5 -i*n I ) ^ _ e (wM-^)^(^))(y-yiv-iww-)^ 

x e ^{u[){Y w -\SY w ) e u>{y' 2 ){Y bl ~\SY bl ) ^ _ e (u;(i/) - u(v[)-u{v' 2 ))(Y bb -\5Y bb j^ (3.2.2.9) 

In the limit that iv ± 1/2 — > 0, then w(i^) — ► —as/(iv ± 1/2) , and in our kinematic region, 
Eq. ( |3.2.2.S| ) reduces to the following expression 



? Y rY b - b -\8Y b 1 f°° ^ /^N** 



dY 2 I v l v J v J'l[ I ' 1 



I Yw+¥ y; L J» Ul U2 Vl V2 " «| ^^±1/2)^^1/2)3 Vfc 2 , 

x r v x L-^d Vl u 2 L-pdu 2 (y-Vw-Ww) {Y-Y bb -±SY bl) 

J— 00 J~ 00 

x U _ ^M-uM-uteMY-Yw-lSYw)^ 

x f°° v[L-/-du[ r v' 2 Ljhv' 2 e^)Va-¥y b i) (Y bb - ¥ Y bh ) 
J —00 1 J —00 2 

x - e( u ^KJ-^KDftr j (3.2.2.10) 



In Eq. ( |3.2.2.10|) , the integration over v can be taken analytically, since 
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1 



(3.2.2.11) 



k^=k 2 



-tYy 



oo-i( e -l/2) du 1 

-oo-i(e-l/2) 2^i (i!/±l/2) (1/2 T^) 3 

oo-i(e-l/2) ^ e ~ i/2±i„ Yv 

-oo-i(e-l/2) 2^7 (^±1/2) (1/2 T ^) 3 

« e((y-yi + y 2 ) + (3.2.2.12) 

+ terms suppressed by powers of as relative to the first term (see Eq. ( 3.2.2.11] )) (3.2.2.13) 

3 

(3.2.2.14) 



1 

a s 



dt 
2nt 




sp 



- exp [2Ja s Y v \n{V/kl)) 1 



a s 



t 



sp 



where t = —as/ (1/2 ± if), and Q(x) = 1 if x > while 0(x) = for x < 0. The last integral is 
taken using the steepest decent method, and t S p = 



a,S ln(fc 2 /fcg) 



Yi 



is the saddle point of the integral in 
Eq. ( 13.2.2.141) . fg P = t S pYi. However, this method works only in the region Y{tsp > 1. For Yitsp < 1, 
we can use Eq. ( |3.2. 2.13| ). Using Eq. (|3.2.2.13|) for the calculation, we find that the largest contribution 
stems from the first term in both parentheses in the integrand of Eq. ( 13.2.2.10 ). Hence, the first terms in 
both parentheses survive, leading to the following expression 



Y 

dYi 

Y W + ±6Y W Jo 



Y bb 5^M> 



dY 2 I u Iu! Iv 2 I' U1 K 2 



CO poo 

vtL-^dux / u 2 L-^dv 2 I u[L r i du' 1 I i/ 2 Ly*di/ 2 X 



oo poo 

,1 r-^J../ / ../ t~ %v 'i,U.I 



x J_ e u(vi){Y-Y w -±8Y w ) e u(u 2 )(Y~Y bi ~\5Y b - b ) ^(^J {Y w -\SY W ) ^(u' 2 ) (Y bb -\5Y bb ) ( 3 . 2 .2.15) 

ai 



The integrals over z/j and z/-, lead to the Pomeron amplitudes given by Eq. (|1|), and the final answer 
for this kinematic region is given by the following expression. 



A (IP ^2P - JP;Fig.$) ieg n 



(27T) 



16 



-E Ik*, M w ,M bb , ho) uf p 4 P vf p v'i p x 



a- 



(3.2.2.16) 



P (k, M b - b \Y - Y bl - ±6Y bb A P (k , M b - b \Y bb - Uy^ P (k, M W \Y -Y w - ±8Ywj P (h>, M W \Y W - ±5Y W \ 



where vf is the value of the saddle point in the integration over z/j, in the steepest decent method, 



and 



..SP SP 'SP iSP 
u l v 2 v l u 2 



7T M M w/k ) HM^/H) imti/n H M b y k o) 

(Y - Y bb - \5Y bb ) (Y-Y W - \8Y W ) (Y bb - ±5Y bb ) (Y w - \8Y W ) 



- 20 - 



3.2.3 Region (III): \iv + 1/2| « v x + i/ 2 + 1/2| « i/( + u' 2 ) 

As it has been pointed out in Ref.[||], the fastest increase with energy stems from the third region of 
integration over v, v\ and v 2l which is specified by 



v 2 ; 1/2 + iv — iv\ - iv 2 



'2 I 



1/2 + iv' 



IV ■> — IV n 



0: 
0: 



u(2v% + \ i) 
u{2vf + U 



2uj{v[°) 



(3.2.3.1) 



It turns out, that 2u>(vf) > 2u;(0), and therefore the contribution from the kinematic region specified 
by Eq. (|3.2.3.1 ), leads to a faster growth than in the case of the exchange of two non interacting Pomerons. 
It means, that this region can lead to a larger value for the cross section for the two parton shower 
production at the LHC energy, in comparison with the Tevatron energy, (where such a cross section has 
been measured in ref. Q). We can integrate over Yi and Y 2 in Eq. flOg ), to give the result 



9 

/oo _ roo poo /'CO 

dv v 2 Yl / dvivf A (vi) / dv[vf \{y[) / dv' v' 2 (3.2.3.2) 

-oo J —oo J — oo J —oo 

x g P (k,q = 0;-v) / d 2 </ r 3JP (g = 0, q '\v, v x , v 2 ) T 3]P (q = 0,q'\ - v' , -v[, -v' 2 ) g P (k , q = 0; v' 



-l 



x o- (Mm, o- (M 6 5, -i/ 2 , ^4) (^(^l) + w(z^) - lo(v)) - lo(v[) - u(v' 2 )) 

x ( e u{vi) (Y-Y W -\SY W ) e u,(v2) (Y-Y bb -±5Y bb ) _ e u{v) {Y-Y W -±5Y W ) e (u{v 2 )(Y w + \5Y w -Y bb -\5Y bb ) 



x y e ^(K)(Yw-±5Y w -Y bb + l6Y bb ) e u>(v>) {Y bb -\5Y bb ) _ e ( u (^) (lV-§*y w ) e «(^) 
Using Eq. ( 3.2.7|) and Eq. Q3.2.9 ), and the expression given in the appendix (see Eq. (|2.12j )), we obtain 



o 

/oo /*oo J' roo roo 

dv / dv' Yl / v 2 A (v^ dvi / v' 2 X (v'i) dv[ 

-OO J —OO J —oo J —oo 



x (16k 2 ) 



'16k 2 



'16M 



2 xJ^-Ji/i 



16M 



66 J 



<5(l/ — V\ — V 2 — v' + v[ + i/ 2 ] 



gwfc) (y-yvF-i«V) ^2) (y-ns-|«iB) - e ^ u) ( Y ' Y w-h 5Y w) (Y w +i5Y w -Y bb -iSY bb ) 
(1/2 — iv + ivi + ii/ 2 ) 3 (w(^i) + lo(v 2 ) - lo(v)) 

(1/2 + iv' - iv[ - iv' 2 f(uj(v') - u(v[) - uj{v' 2 )) 
xdi-vj) d(<) C x {-v 2 ) C x (v' 2 ) 



(3.2.3.3) 



In order to evaluate the integrations over v\, v 2 , v[ and v 2 on the RHS of Eq. fl3.2.3.3Q , it will be 
useful to make the following change of variables. Let v>\2) = 5 ( u l + ^2) > an d ^[12] = 5 (^1 — ^2) , and 
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likewise for v[ 12 \ = | ( v\ + v%) , and z/ 12 i = | ( i/[ — i/' 2 ) . In new variables Yii=i v l ^ (^) ^ /2 ^ W) ^ 



^(12) ^[12] ^(12) ^'[lz]' ^ e ex P ec ^ that V{Y1) anci ^'i^ wm 



(12) 



'(12) - 2 VI T- "2; , ^ [12] 

has the form 16 (v 2 (l2) - i/p^) (i/g 2) - z^[ 2 12] 
be close to z/^ while both z/j 12 ] and z/ 12 i will be much smaller. Indeed, assuming this we can expand 
uj{vi) = w(V(i2)) + l ^ / ( l/ (i2)) ^[12] ~~ w "( l/ (i2)) v \<&l^" Using this expansion one can see that integrals over 
v\\2\ an d ^[12] have a Gaussian form and can be taken using the steepest decent method. The saddle point 
value of these variables are small and decreasing with energy, namely, 



SP 

[12] 



In {M&m/k® 



2Y - Y w -\SY W -Y h 



bb 



\SY t 



JSP 

'U21 



In (M& Ml/kt) 



bb 



Y w - \8Y W + Y b 



bb 



\8Y> 



bb 



(3.2.3.4) 



Evaluating the integrals over V(y2) an d v '(i2) ^ closing the contour on the pole, and by using the 



8 function, Eq. ( 3.2.3.3 ) reduces to the following expression (considering i/r^i as well as z/^j are small: 



.SP 



V [12]) < V {V£) and ^[ S 12]) < ^(12))- 



A(jT -> 2W> -> JP; §) = vr 8 £ (fe, M w ,M bi , k ) J dv dv' v x v 2 dv 12 v[v' 2 dv' l2 



M w M l, 



{uj(vi) + oj(v 2 ) - w{v)) (uj(v[) + w(f 2 ) ~ 

x ^(i^) (y-y^-inw) gwc^) (y-yrt-s^k) (i - e M^-"(n)-^ 2 ))(y-yw-f5yw)^ 

x ^(^(Yw-lSYw) e u>{v' 2 )(Y bl -\5Y bl ) f x _ e («(i/)-a>(i/ 1 )-w(«4))0kk-3* y tf) > 



(3.2.3.5) 



In this integral zv. 



(12) 



I ») and "(12) 



\W ~\ *) an£ i "i = \(v — \ i) ± v 



121 



i/? at 



high energy, as well as = g(z/ — | i) db z/ 12 i — > ^ . It should be stressed, that in Eq. ( |3.2.3.5| ), all the 
z/s satisfy Eq. ( |3,2.3,1| ). Expanding uj{vj) around the point in the vicinity of v\ = u[ = v®, and denoting 



iui — %v\ = f and iv\ — iv\ = /', we find that Eq. (3.2.10) can be written as follows 



e — ioo 
e+ioo 



(u(^) + d//2) (2Y - Y w -Y b - b -\(SY w +SY b - b )) _}_ ( j _ d / (r-y^-i^A 
2tH d/ V / 



( w K) + d/'/2)(y w +y i , 5 -i(5y M +<5y i)5 ) J_ , , 
Tan df 



-df (Y bl 



(3.2.3.6) 



where d = 2a;" (v\ 



> and d = d — ll>" (2v® + il / 2 > d. In this integral, we can close the contour 



of integration over / and /', on the poles / = and /' = 0, for the first term in each set of parentheses. 



In all other terms in Eq. ( 3.2.3.6 ), we can close the contour on the right semi-plane, where we have no 
singularities, and therefore these integrals are equal to zero. Finally, we have 
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A(JP 2iP IP; Fig. 

K) 4 r dv^dv\. 



d 2 



'12 a^n 



Tr s Z{k,M w ,M bl ,k 
k 2 



(3.2.3.7) 



2u° 



7.2 




2v? 



w"(o) (y - y w - 



xe 



( in Eq. ( |3.2.3.7[ ), we hope that the factor d will be not confused with the sign of the differential) 
Expanding oj(y^ + \vvi) and w(y\ + ^u[ 2 ) for the case of small (v'12), and integrating over 7^12 and v' 12 , 
we obtain 



A (IP ^ 2iP -» iP; Fie/. 



7r 9 S(A;,M w ,M 6 6,fco) 



0\4 



K) 
d 2 



x w "(0) (F - y w - ijYiy) exp (2 u(yl) (Y - l -(5Y w + SY bb )) 



^2lu"(0)(2Y -Y w -Y b 



exp 



In 2 (M w M bb /k 2 



bb) 



2u>»(0)(2Y-(Y w + Y a )) 



V2w"(o)(y w + y t 



: exp 



bb) 



In 2 (M^M^/fc 2 ) 
'2a;'' (0) (Y w + Y bb ) 



k 2 



2v° 



1- 





2v° 



(3.2.3.* 



Eq. (3.2.3.8) can be rewritten in a simpler way, namely 



A (JP -» 2iP iP; D m = vr 9 E (A;, M w , M fe5 , k 



x o/'(0) (y - Y w - ^SY W ) exp ( 2 



,0\4 



d 2 



2/' V 



/• 2 



2^? 



^(0)) (y - \{5Y W + 6Y bb ))\ \ (3.2.3.9) 
x P (k, y/M w M bb \2Y -Y W - Y bb - ^(5Y W - 5Y bb )^j P (k , ^M w M bb \Y w + Y bb - ^(5Y W - SY bb )\ 

3.2.4 Region (IV): u < 1, v x < 1 and v 2 < 1 (1/ < 1, v[ < 1 and v' 2 < 1 ) 

In our previous calculations, we focused our efforts on the singular part of r^p (q = 0, q'\v, u\, V2) = 



2 ™s 



J oP/c g(k,q = 0, v) g [k + ^q',q',vi) g ( & + \q' ■, — q', v 2 I (see Eq. (|3.2.8D ). However, the inte- 



gral in the region where all the t/s are small, also can lead to a contribution which increases at large values 
of y, as exp (tu(0) (2y — ^5Yw — \5Y b if). Indeed, in this kinematic region we can integrate Eq. ( 3.2.1.2 ) 
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over all z/j (v[), using the expansion of Eq. ( |2.15| ), As an example, the v\ integration in Eq. (3.2.14) takes 
the following form, using the expansion of Eq. Q2.15 ). 



dv! vtX (1/1) e ^(o)(n-yw--|5y w >i 1 , 1 V'(o)(y 1 -y w --|5iV)^iin(Af2,/fc 2 ) 



(3.2.4.1) 



The exponential has a saddle point at the point z^ p = iln {M w /k"^\ / (u" (0) (Y\ — Yyy — ^SYw))- 
The expression in the exponential / (t'l), at the saddle point takes the value 

f(uf) = ±ln 2 (M^/fc 2 ) / (u"(0)(Y 1 -Y w -lSY w )). Hence, evaluating the integration in Eq. ( |3.2.4.1| ), 
using the method of steepest descents gives the following result. 



16 ^(0)(^-^) x i-^]' 2 = 16- f ^ 

V /" K ) / ^" (°) ( y i ~ Y w-¥ Y w)) 



< e u{0)(Yt-Y w -\5Y w ) x e l\n 2 (Ml/kf ) )(u J "(0)(Y 1 -Y w -^SY w )) 1 



2tt 



u"(0) (Yt-Yw-lSYw) 




where /" (u[ p ) = uj" (0) (Yi -Y w - \5Y W ). The factor of 16 in front in Eq. ( p.2.4.2j) , originates from 
the value which A (^i) appearing in Eq. ( |3.2.4.1| ) takes at small v\, (see Eq. ( |2.8| ) ). The integration over 
the v\, v[ and u' 2 , are evaluated in a similar way, using the method of steepest descents. Hence, the result 
after integration of the RHS of Eq. ( |3.2.1.2j ), can be written as follows 



A (IP -» 2P -> P- Fig. |) reg (iv) = (I6vr 4 ) S(A;,M w ,M b - b ,A;o) (3.2.4.3) 

„y „y-_I,5y- 



2^ X 1 / 2 / 2^ ^/ 2 



In 2 (M w /k 2 ) In 2 (M 2 5 /fc 2 ) 



x 



' \cu"(0)(Y -Y 1 )J \u»(0){Y 2 -0), 
2vr L"(0) {Yx -Y w - ^6Y w ^j u" (0) \Y W - Uy w - Y^j \ 

2vr L"(0) (y x - Y bb - l -5Y b ^j u" (0) (V 6 - 6 - - Y 2 



In 2 (M^/fc 2 ) \ / 1 In 2 (M^/fcg) 

6XP ^ 2 w » (0) (Fx - *V - \ 8 Y w ) J exp y 2 u» (0) (Y w - \ 5 Y w - Y 2 ) 

x exp — — — = r- exp 1 



" 2 u» (0) (Yi -Y bb -±5Y bb ) I 2 u/' (0) - 1 5 Y fe5 - Y 2 ) 
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The largest contribution to the integrals over Y\ and Y2, stems from the region Y — Y\ oc 1/^(0), and 
Y2 — oc l/ui(Q). Since a;"(0)/o>(0) « 1, as far as the QCD coupling is concerned, we can conclude from 
Eq. ("3.2.4.31) that (using the Eq. (2.17) for the expression for P (k, fco| Y — Y'), for the case of small v ) 



A(P ^2P -> P; Eq. ( ^2~4~3l) ) rog T v oc (2vr) 16 S (fc, M w , M bb , ko) (3.2.4.4) 

vr 2 ln(M^/fcg) ln(M%,/fcg) ln^/fc 2 ) ln(M£/fcg) / 1 

X (y " % " W*) (Y-Y W - \5Y W ) (Y bl - \5Y bl ) (Y w - \SY W ) \ ' b > bl bb Y bb 

x P (ko, M b - b \Y bl - Uy^ P (k, M W \Y -Y w - \sYwj P [ko, M W \Y W - ±6Y W ) 

Comparing this contribution with Eq. ( |3.2.2.16| ), one can see that 

A(P^2P^P; Eq. (jUg )) jy 
A(P^2P^P;Eg. (|3.2.2.16|) TT ^ ^ « 1 

Therefore, this kinematic region gives a much smaller contribution, than the kinematic region (II). 
4. Estimates of the background 



Comparing Eq. ( 3.2.2.16} ) with Eq. ( 3.2.3.9| ), one can see that the first contribution is suppressed. Indeed, 
the ratio of these two contributions for Yw = Y bb = in the c.m.f., is equal to 

A (jP -> 2P -> P; Eq. Hp reg III _ 2 2 VwWl72 1 5Y ) K) ' 

A(P^2P^P;E q . (b.2.2.ld) )vezTl ~ ^ ( ° } {Y/2 " 2^ (F/2 " 2^ * v sp y SP u[ SP y ,SP 

x exp (2 [u{y\) - u(0)) (y - l -(5Y w + 5Y a ))j ) » 1 (4.1) 

Indeed, the main factor that determines the energy behavior, is 
exp (2 — w(0)) (Y — h(SYw + SY bb )) ) 1 which increases with energies. Other factors include 

cj" 2 (0)/d 2 w 1/4, and a 5 (F/2 - §<^V) > 1, as well as a s (Y/2 - \5Y bb ) > 1. However , at high energy 
v\jvf p S> 1. This factor cannot change the asymptotic behavior of this ratio, but it leads to a decrease 
in the range of energy that we are interested in, namely from the Tevatron to the LHC energy. 

It turns out, that both at the Tevatron and the LHC energies, the ratio of Eq. ( [O]) is smaller than 
unity , since vf p is equal to 0.37 and 0.19, respectively, if we take as = as (Mh ko), and k$ = IGeV. 
Such values of i>sp lead to v\jvf p = 0.21 at the Tevatron, and v\jvf p = 0.42 at the LHC. The numerical 
solution of Eq. (p.2.3.1| ), gives the following relation (see Ref. ||) 

2{u{y\) -u>(0)) ~ 0.25 a s (4.2) 
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which leads to exp (2 - w(0)) (Y - \{8Y W + SY$))) ) « 1.4(for the Tevatron and 1.86 

for the LHC. Collecting these, estimate we obtain that the ratio of Eq. (|4.1|) to be equal to 7 x 10~ 4 for 



the Tevatron, and 1.4 x 1CT 2 for the LHC. It means that we can neglect the contribution of Eq. ( |3.2.3.9|) , 
and consider only the contribution given by Eq. (|3.2.2.16|) , which describes the production of the W boson, 
and the 66-pair from two parton showers. 

In this paper, we obtain two different expression for the two parton shower production, namely 
Eq. (|TD for the non-enhanced diagram, and Eq. ( |3.2.2.16l) for the enhanced diagram. From a purely 



theoretical point of view, it has been proven []T3|, that the only contribution stems from the enhanced 
diagram, which reproduces the non-enhanced diagram in the low energy limit (Y — y\ ~ l/as,V2 ~ ^-/ a s)i 
at least for the case of dipole-dipole scattering. In this case, their contribution given by Eq. ( p. 2. 2. 16 ), has 



been measured at the Tevatron ||, and it's value at the LHC energy can be estimated using the following 
expression 

d 2 a + p -» W + [66] + X) 1 uf p (LHC) i/f p (LHC) v[ SP (LHC) u' 2 SP (LHC) 



dYw dY bb ^ a eff vf p (Tevatron) v 2 p (Tevatron) u'f p (Tevatron) v' 2 SP (Tevatron) 

r'2 



da md [ln( S /M b - b y,Y bl ;0) d (ln^/M 2 ,); Y w ; 0) 
X dY^ dYw (4 ' 3) 



! da ind [Hs/M^); Y bl ; J d(Tmd (hi^/M 2 ,); Y w ; 0) 



16 a e ff dY bb dY w 

The ratio of this background to the signal, can be written as follows 



M, 



dY w dY bb __ dYw 1 S/f 2 d a incl (p + P^bb + X) 



M - v - >- (4 4) 

d?ir(p+p->w+H+x) <P<r incl (p+p-+w+H+x) iQa pff bb dM 2 - dY h - h y ' ; 

dY w dY H dY w dY H ejJ bb bb 

The first factor in Eq. ( ^4.4| ) is well known (see Eq. (|1,1|)), and the uncertainties in the values of the 

gluon structure function does not affect it, to within an accuracy of 7%. To evaluate the cross section 
M?r d — i we rewrite this expression as follows 

bb dM[\ dY bb ' * 



2 d 2 a ind (p + p^bb + X) _ a (M 2 g ) da (p + p -> H + X) 
bi dMl dY bb " a (M 2 ) dY H [ " 

The values of a (M 2 ^) and a (M|) are calculated in Appendix 2 and Appendix 3 and it is equal to 

rr ( M 2 -) A QO. 

s m(M 2 5 /4m 2 ) « 357 (4.6) 



°- «) 4 32a 2 

^ 3M^ " 
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We believe that this ratio, which we obtain in LO QCD, does not depend on the accuracy of our 
calculation, and will be almost the same as in the NLO approximation of perturbative QCD. All the 
uncertainties related to the calculations in high order perturbative QCD in Eq. fl4.5|) , are absorbed in the 
value of the inclusive Higgs cross section. For this value, we take 3 x 10 -9 mb (for LO QCD), and 5 
x lCT 9 mb (for NLO QCD) (see Refs. Substituting this value and Eq. (fl|) into Eq. (U), we obtain 



r, d 2 (Ti nc l (V + P — > bb + X) n , 

< IM^ 1 = + 2 > * 10 ~ 8 mb (4 ' 7 > 

Using this value, we estimate the background-to-signal ratio, as 

M bl dY wd y b - b 1.24-2.4 



da(p+p->W+H+X) IQ 
dY w dY H 



Therefore, if Eq. ( |3.2.2.16 ) is responsible for the two parton shower production, the background in the 



LHC kinematic region will be small, and it can be neglected. 

However, if the proton has some non-perturbative contribution, and differs from the colourless dipole 
of perturbative QCD, we can hope that the non-enhanced diagram gives the main contribution, leading to 



M, 2 7 



d 3 a{p+p^W+[bb]+X) 



dY w dY H 

The background of Eq. (|Q|), should be considered . However, using the fact that the production 
from the two parton showers does not depend on rapidity, while the signal leads to the sharp rapidity 
distribution in the rapidity of the bb pair, the background can be easily separated. We would like to draw 
your attention to the fact that all our estimates, are heavily based on the beautiful measurement, of the 
so called double parton cross section, by the CDF collaboration. We would like to draw the attention of 
experimentalists, to the fact that this kind of measurements in the LHC region, will clarify the value of 
the contribution of the two parton shower processes, and therefore, will lead to an improvement of our 
knowledge of the non-traditional QCD background. 



5. Summary 

This paper has two main results. The first, are the estimates of the background to the associate W Higgs 
production. The second result, is the detailed calculation of the contribution of the enhanced diagrams 



to this background. As far as practical estimates are concerned, Eq. (;L8) and Eq. (L9) give the main 
conclusions of this paper. These results are encouraging, since the background turns out to be no larger, 
than a factor 2 larger than the signal. The principle difference between the signal and background, is the 
dependence with respect to the rapidity difference Yyy — Y^. The signal has a peak at Yy/ — Y b i = 0, 
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while the background does not depend on this difference. Having this in mind, we can easily separate this 
background for the signal. 

The theoretical calculations, are the first where all kinematic regions in the enhanced diagram have been 
considered. The most interesting part, is related to the calculation of the contribution of the overlapping 
singularities. Being interesting theoretically, these singularities give a negligible contribution at the LHC 
energies and below. We believe, that this is an important observation for practical estimates of the 
processes, that are responsible for physics in the saturation domain. 

We would like to draw the attention of experimentalists, to the fact that the measurement of the so 
called double parton cross section in the LHC range of energies, in the way that the CDF collaboration 
did at the Tevatron, will clarify both the contribution to the background of such non-traditional QCD 
processes as the two parton shower production, and the importance of the shadowing corrections for 
semihard processes. 

We hope that our paper will be useful for planning experiments, in the search for the Higgs boson. 
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A. Appendix 

A.l The amplitude for the subprocess in GG — > W 

The aim of this section of the appendix, is to derive the expression for the cross section a(M b i), for the 
production of the bb quark - anti quark pair, in the gluon fusion GG — > bb. These cross sections, have been 



calculated (see Refs.[17, [j^, 19 1), but we reproduce the calculation here in our technique using Sudakov 
variables, since using this method of calculation, it is easier to include the calculated cross sections in the 
framework of our approach, based on the BFKL Pomeron. 

In Fig. [|, one can see that in the Born diagram for this process, the hexagon is on shell, as well 
as the protons between the t-channel gluons. In other words, the dashed line through the center of the 
diagram, intersects lines representing particles which are on shell. The cross section which corresponds to 
this diagram, is equal to 
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da BA 1 (\ \ 

— = — A (Mueller diagram of Fig. |5| for process p p W — ► p p W) 8 I - In (ai//^) — yw 

aYw 4s 1-1 \2 / 

(A.l.l) 

where a and (3 are the Sudakov variables for the vectors q% and q\ ( see below). The upper line of the 
diagram gives the factor 

(-^r) u(pi, Sl )^u(pi -q x ,s[) x (^p) u{ Pl - gi , s[) Yu{ Pl , Sl ") (A.1.2) 

where si, s' x and si" are the helicities of the incoming, the intermediate, and the outgoing quarks 
respectively. Throughout this calculation, the Eikonal approximation is used, which assumes that the 
components of the momentum of the exchanged gauge particle, (in this case gluons), are small compared 
with the momentum of the emitting quark. Hence, in the Eikonal approach, all components of q\ are 
assumed to be small, such that Eq. ( |A.1.2 ) may be replaced by 



/ 47rag \ 
V 4 J 



(pi, si) Y l u( Pl , s[) x u( Pl , s[) -f p u(pi , si") (A.1.3) 



For spinors which are normalised such that u (p± , s^) u> (p± , si) = ^s^s^, and using the following 
property of Dirac spinors from Gordon's relation, namely 

2u(p, s) 7^n (p , s') = u(p, s) ( (p + p'Y + io^ (p-p') u ) u(p, s) (A. 1.4) 



then plugging this into Eq. (A.1.3), one has for the contribution of the upper line of the diagram 



(Anas) M o x x (4vra 5 ) Qi,± , . 
i — PlPl 5 s u s' o s ', Sl " = 4 o > 5 S > n (A.1.5) 

where we use the gauge invariant condition, that q± H^^^^ = 0, where H is the amplitude of 
the hexagon. Similarly, the lower line of the diagram in the Eikonal approximation, gives the following 
contribution to the amplitude 

(4vra s ) v a , x (47ra 5 ) q% t ± Q° 



Hence, the expression for the amplitude of the diagram shown in Fig. [5], is given by the following 
expression, where it is assumed that the on shell quarks (horizontal lines in the hexagon shown in Fig. ||), 
carry momenta which depend mostly on the longitudinal direction, and all other quark lines in the hexagon, 
carry momenta which depend mostly on the transverse direction, 
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47rag (g 2 )) 2 t d 4 k f dH 
\h\oc\j5\ / (2tt) 4 7(2^ 



^(*%-§ = V 4 4 2*2 /7^4 /t^T# ( A - L7 ) 



x (27r) 2 6 y - (3 k ai s- (%-qiJ ^- m 2 ^ 5 ^ a//3 2 s - (l + <f 2 ) ^ - m 2 d 



where H = ± , q[ ± , q% ± , q% ± H j 



fipua 



1 { 1 hv) - e Vtr ( 7A ( 1 - 75 ) (/ - ^) g^ 7r ( 1 - 75 ) fy 2 (> + 



fl g 2 gfe Vtr ( ^//^ (1-^) *) 



'(%) o 2 2,2,2 
"1P2S 9l ^ ' 



and where cr ( M w ) = 4 ( 4vra s (M w )) 2 V2G F 1 ' 



N 2 -l 



(A.li 

2 -, 



2A^ C 7 4iV c 
Z?i = A; 2 — to 2 = OLkflk s — k\ — m 2 

D 2 = l 2 -m 2 d = aifii s - li - m 2 d (A.1.9) 

where in the last step, 7 5 anticommutes past 4 gamma matrices in the trace, to give a ( 1 — 7 s) 2 = 
2(1 — 75) term. We use the fact that the part of the trace on the RHS of Eq. ( |A.1.7| ), which contains a 
7 5 matrix can be removed, because evaluating the part of the trace containing the 7 5 matrix gives a result 
which is antisymmetric in fi and p, which contracted with the symmetric term p^p^, gives zero. e A and e A 
are the polarisation vectors of the final state W bosons in Fig. ||[ In Eq. (|A,1.7| ), we use the fact that 

eV = g Xr (A.1.10) 



and using the integration measure 

d A k f dH sT. .„ d 2 ki s r , .„ d 2 l 



Jda k d(3 k ^ S - (A.l.ll) 



(2vr) 4 J (2tt) 4 2 J rK (2vr) 4 2 7 "* (2vr) 4 

then Eq. ([A.1.7 ) simplifies to 



A (Fig. 



(47r«g (g 2 )) 2 
x 5 I — /SfcCKi s 



ci 2 fc 

(27T) 



L £ 
3 2 



/ 



dai dpi 



d 2 l ± 
(2tt) S 



<5 a ; /3 2 s - ( I + g 2 



(A.1.12) 



- 30 - 



Using all the 5 functions to evaluate the relevant integrations, specifically taking the integral over 
ot\ > ai and over (3 2 > Pk, and using Eq. ([A.l.ip , we obtain 



da BA al(q 2 ) x 

—— = 2 syq 2 > a(M w ) 

dY w ii±qi± 

It should be stressed, that the factor a 2 s (q 2 )/ (q 2 j_ q\ jY is included in the exchange of the BFKL 
Pomeron, or in more general words, it should be absorbed in the product of the two gluon functions 
x\G{q 2 , x\) X2G(q 2 , X2) that give the flux, (luminosity) of the gluons from both protons. The origin of the 
logs in Eq. ( |A.1.13| ), comes from the integrations over k± and l±. These integrations take the following 
form 



/ 



d 2 k±k 2 ± ( from the tr(. . . )) — = vr ln(M^ / '4ml) (A.1.14) 



A. 2 The amplitude for subprocess GG — > bb 

In this appendix, we calculate the amplitude for the subprocess GG — > bb, in the region of large values 
of 66 mass, ( denoted by M h i « M#). For such large values of the 66 mass, we can restrict ourselves by 
finding the contribution of the order of In (M^/4m^) , which stems from the simple diagram of Fig. [6|. The 
cross section can be calculated using the following expression 



d 2 a BA 1 



M^g 2 = - — A (Mueller amplitude (see Fig. |6|) for the processes p + p + [66] — > p + p + [66]) 

"■ M bb° bb 4s 

x Ml 5 ( ai (h s - Ml) 6 Q In - Y bI ) (A.2.1) 

The amplitude for this diagram is given by the expression 



A (Fig.® 



d^k ( ig s ) _ . . „ / , r\ ( Ws) - 1 \ v ( „/ / \ /aoo\ 
~i — u{Pl, si) Tu( Pl Sl ) —u{p 2 , s 2 ) 7 u(p 2 s 2 ) (A. 2.2) 



(2vr) 4 q{ 



<l-2 



u{Pi, si) Yu (p' x s[) y % 2 s ' u{p2 , s 2 ) l a u(p 2 s 2 ) ) S", 



9? 



t 



Using Eq. ( [AToD and Eq. ( jAXg ), we find 
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A(Fig.§ 
where S 



4iras 47ras f d 4 k 
(A^a s (M bl )) 2 



S 



(k 2 



m 



^i,±^i,^ g 2,±92,± tr { (ft - i x + m b) 7 P (ft + m b) la (ft + ^ 2 + rrib) 7u (ft + m 6 ) 7 M | 



x (2vr) 2 J^fc-gi) 2 -m^ <J HA? + g 2 ) 2 - mft) 
In Eq. ([AlDD , we have used the Sudakov expansion 



N 2 -l 
N c 



V 



9; 



k+q. 



Figure 9: The process G + G — ► 6 + 6. 
following expression 



<?2 



aiPi + P1P2 + qu 
U2P1 + P2P2 + <?2J 



(A.2.3) 



(A.2.4) 



To calculate the trace in the expression for S given in Eq. ( A.2.3 ), 
we recall that at high energy, the horizontal on mass shell b quarks 
in Fig. ||, have momentum which depend mostly on the longi- 
tudinal direction, and the vertical virtual b quark propagators, 
have momentum which depend mostly on the transverse direc- 



tion. This observation reduces the RHS of Eq. (A.2.3) to the 



tr I (ft - f x + rrifc) $ 1± (ft + m b ) ^ ± (ft + f 



+ q 2 + m b 



$ 2L (ft + mb)$ ls \ = Aa 1 f5 2 sqlql(k 2 +ml) 

(A.2.5) 



In terms of the Sudakov parametres introduced in Eq. ( A.2.4j) , the integration measure of the k 
momentum space, becomes 



J (2vr 



4 k s f d?k±dpd\ 
(2ttY = 2 J (2tt) 4 



(A.2.6) 



Using Eq. ( |A.2.6| ), we integrate Eq. (|A.2.2Q over all the as and /3-s, taking into account the 5-functions 
in Eq. ( A.2.1|) . The result is given by the following expression 



„ d 2 a BA a 2 (q 2 ) 

bb rfTU-2 JV_ „2 2 V obJ 

aM bb a bb Ql,± %,±_ 

where „ (A % ) = (%T) ' -J- 5* In (>% ) (A.2.7) 
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where the log term has the same origin as in Eq. ( |A.1,14| ). 



A. 3 The amplitude for the inclusive Higgs production via subprocess GG — ► H 



The expression for the cross section for the inclusive Higgs production process shown in Fig. 10, is very 
similar to the process of the W production due to gluon fusion, and it is equal to 



da 



BA 



— A (Mueller diagram of Fig. |To| for the process ppH 



p p H) 5 ( -ln(ai//3 2 ) - Vh 



(A.3.1) 



The amplitude for this process shown in Fig. 10, is given by 



A (Fig. M 



a n^jyiipo- 



(47ra 5 (q 2 )Y 

4 4 



(A.3.2) 



x ^ {u (p[) 7 M n (pi) u (p' 2 ) Yu (p2)} {u (pi) Yu (pi) u (p' 2 ) 7^ (p 2 )} 
spins 



where A 



AV2G F a 2 s (M H ) 



1 



and 



v M 
QlQ2 



(A.3.3) 





Figure 10: The process of inclusive Higgs 
production due to gluon fusion: G + G — > 
Higgs. 



Using Eq. ( |A.1.5| ) and Eq. (A. 1.6), as well as the form of 
the gauge invariant tensor H^ u , we readily obtain the following 
expression after integration over the as and /3-s 



da 



BA 



where 



dY H 
a(M H ) 



N?-l\ A 2 



4N 2 J 256vr 



(A.3.4) 



Since the Pomeron exchange, or more specifically the gluon 
structure functions enter in the same way for the case of the 
bb pair production, as in the case for Higgs production due to 
gluon fusion, we can calculate the ratio 



2 d 2 a(p+p- 



-x) 



bb dMfadYn 



a(M l 



bb) 



a(M H ) 



(A.3.5) 
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Using Eq. ( |A.2.7| ) and Eq. ( [A.1.13P , we have 



A [p+ lbi]+x) = 4 *r<p + p~H + X) 

where the factor of 4/3 is a colour factor. 
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